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Abstract 

We study field configurations in a hot quark-gluon plasma with spherical 
symmetry. We show that the electric fields point into radial direction and solve 
the effective non-abelian equations of motions. The corresponding charge 
density has a localized contribution which has a gauge invariant interpretation 
as a pointlike color charge. We discuss configurations oscillating periodically 
in time. Furthermore, we calculate the electric field induced by a constant 
local charge that is removed from the plasma for t > as a model of a decaying 
heavy quark. 

ll.lO.Wx, 12.38.Mh 



Typeset using REVT^ 

■te-mail: herby Iapphpl.in2p3.fr 

*URA 14-36 du CNRS, associee a I'E.N.S. de Lyon, et au L.A.P.P. (IN2P3-CNRS) d'Annecy-le- 
Vieux 

1 



INTRODUCTION 



Unlike as QED where linear-response theory supplies a straight-forward definition of the 
Debye-screening potential as the energy induced by a static delta-like charge distribution, 
the non-abelian analog is not a gauge invariant way to couple a classical charge to the gauge 
potentials. Making the perturbed Hamiltonian gauge invariant generally amounts to include 
ghost fields which compensate for the trivial transformation behavior of the classical field 
strength. It was argued that particularly in the temporal gauge {Aq = 0) the gauge 
potentials and fields are related in a quasi-abelian way hy E = —d'^A, rendering the abelian 
response-theory applicable. 

The TAG, however, does not admit an imaginary time formalism (ITF) in the usual way 
0, and a reformulation in the so-called physical phase space formulation ^ again would 
call for the inclusion of vertices in the linear response function necessary to account for the 
non-abelian character of the theory. 

In order to avoid these difficulties one may consider the electric dispersion relation 
without external sources and define the screening potential as the Fourier-transform of 
the static propagator on the physical dispersion. I want to point out that taking the 
static limit Il{kQ = 0, k) alone does not define a gauge invariant potential, and thus 
a physical quantity. Although the leading order static self energy is just a constant, 
Ilei{ko,k) = m? = g'^T^{N + Nf /2)/3, a next-to leading order definition requires the in- 
clusion of the physical dispersion. In the case that the propagator had a pole well separated 
from eventually cut contributions, the pole really corresponds to a physical observable [Q. 
However, due to mass-shell singularities appearing within the standard resummation scheme 

the longitudinal self energy turns out to be sensible to a - perturbatively uncalculable - 
magnetic mass 0. Although the electric self-energy does develop a gauge invariant pole at 
next-to-leading order, it is not clear to what observable this quantity corresponds. 

In this paper we study spherically symmetric field configurations at leading order which 
have an interpretation as generated by a local charge in a gauge invariant way. 

I. OBSERVABLES FOR SPHERICAL SYMMETRY 

Let us reconsider the generic idea to put a pointlike static quark in the hot quark-gluon 
plasma. A classical source may be characterized by the following properties. 

(i) The source is assumed to have no internal structure such as spin which would involve 
a magnetic moment. 

(ii) The source is localized and pointlike, and its location is static with respect to the 
rest frame of the plasma. Since at finite temperature essentially infrared physics is 
concerned, one can in fact exclude the location of the charge as singular region from 
space-time. To make this statement more precise, we consider physics only at distances 
greater or equal T away from the hypothetical classical source. 

From (i) and (ii) it follows that observables have to be spherically symmetric. Let us 
illustrate that those two assumptions are sufficient in classical electrodynamics to single out 
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the Coulomb potential as the only solution, up to a multiplicative factor. The fields satisfy 
the sourceless Maxwell equations, 

dE = 0, dxB-E = 0, d x E + B = 0, OB = (1) 

and are spherically symmetric. Due to the famous theorem saying that 'you cannot comb 
a hedgehog', spherically symmetric vector fields can always be written as a gradient of a 
scalar in three dimensions, to wit 

E(r, t) = -E(r, t) E(r, t) = df(r, t) 

r 

B{r,t) = ^B{r,t) B{r,t) = dg{r,t). (2) 

It follows from the Gauss law that 

A/(r,t)=0 f(r,t) = ait) + ^-^ E{r,t) = h{t) d-^ (3) 

and from the second Maxwell-equation m. E = d x B = d x dg = thus h{t) = and 
the electric field E cc f/r^ corresponds to the Coulomb potential. 

Turning to the non-abelian case, the symmetry properties (0) cannot be imposed directly 
on the fields, since those are not gauge invariant quantities. According to their transforma- 
tion behavior we are able to distinguish two classes of quantities in a gauge theory, these 
which transform inhomogeneously as the gauge potentials, 

^ ^ U-'A^U + -U^^d^U (4) 

and those that transform homogeneously 

E U-^EU (5) 

as the electric and magnetic fields, the charge and the spatial current. That quantities will be 
called pre-gauge invariant since the trace of products of them defines observables. Operators 
that have positive and gauge invariant eigenvalues correspond to physical quantities, and 
the eigenvalue is the physical observable. Considering e.g. the eigenvalue problem of the 
Lie-group valued electric fields 

E{x, t) \E{x)) = E{x, t) \E{x)) (6) 

where the hermitian field operator E{x,t) is taken in the Heisenb erg-represent at ion, it is 
obvious that the eigenvalue is real and remains unchanged under unitary (gauge) trans- 
formations (^. In that sense all pre-gauge invariant quantities potentially correspond to 
physical observables. In order to construct a set of equal time observables, however, we 
need those operators to commute with each other. As concerns the Lie algebra su(N) one 
may chose a Cartan-Weyl representation with the basis elements T'^, s = 1, . . . , r, r = 
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rank of the algebra, spanning the abehan Cartan subalgebra and the remaining n — r gen- 
erators T^"" where a is an r-dimensional root vector. A set of commuting variables is given 
by the diagonal components E'^T'^ (no sum over s) and the hermitian commutators 

[E^'^T", (E*'«T")t] = E^'^iE''")* J2 «(^)^' (no sum over i and a) (7) 

s 

where a{s) is some weight depending on the group under consideration. Imposing spherical 
symmetry restricts the real diagonal components E^ analogously to the abelian case (^, 
whereas only the modulus of the complex components E"' of the roots is subject to the 
symmetry, 

— ♦ 

= -|E'^(r,t)|e'^"('-'^'*). (8) 
r 

This can also be seen directly by considering the eigenvalues of the eigenvalue problem (^ 
where only products like E°'E~°' = E°'{E°')* enter. Moreover, on account of the commutator 
[T'^,T"] = a(s)T", the phase of one additional non-diagonal color component in (^ can be 
absorbed with the help of the local gauge transformation U = exp(iT''x"/a(s)) which leaves 
the abelian part i?*T* invariant. This shows, that in a spherically symmetric situation, only 
r + 1 {N for SU(N) ) components of pre-gauge invariant variables can been chosen to depend 
on the radial coordinate and time, but the vector component of the variables has always to 
point in radial direction as in the abelian case. 

One could also consider the more general situation where external space-time and internal 
symmetries mix. For SU(2) as internal symmetry group, this leads to an enlarged number 
of solutions of the Maxwell equations 0, but is in principle restricted to the case where 
internal and external symmetry groups are locally isomorphic i.e. 0(3) ~ SU(2). 



II. SPHERICALLY SYMMETRIC FIELD CONFIGURATIONS IN A HOT 

QUARK-GLUON PLASMA 

Let us now turn to the similar situation in a hot quark-gluon plasma. For field strengths 
of order of the temperature, A ~ T, and in the long wavelength limit, d ~ gT, the Maxwell- 
equations without external source generalize to [§] 

[^.,^'^1=Jind (9) 

where = + ig [A^^, . ] and F^'^ = [D^jD"] /(ig). In that regime the full non-abelian 
character of the theory enters, and one cannot expand the covariant derivative in a series in 
the coupling constant. The induced current j^"^^ describes the response of the plasma to the 
non-abelian gauge fields and has to be a covariantly conserved quantity 

[D„f] = (10) 

for consistency with the Maxwell equations (^. It can be expressed by means of an auxiliary 
field W^' 

jU^)=m' f ^v^W\x,v) (11) 
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where = {l,v) denotes a null vector and the angular integral J dQ runs over all directions 
of V. is a solution to 



D, W^"] = F^^v^ (12) 



and conjugate to v^, v^W^ = 0. The set of equations (y,|ri|) and (|T2D forms a complicated 
highly non-linear system of integro-differential equations describing the long wavelength 
behavior of a hot quark-gluon plasma. A particular class of plain-wave solutions where the 
gauge potential depends on space-time only as a function of {p ■ x) has been found in . 
Here, we will consider spherically symmetric configurations with vanishing magnetic fields 
= e'^^Fjk = 0. Note that under a gauge transformation B transforms homogeneously as 
the electric field (^ so that solutions with vanishing magnetic fields have a gauge invari- 
ant meaning. Moreover, those configurations are singled out in the sense that the spatial 
components of the corresponding gauge potentials can always be written as pure gauge, 

A=-V-^dV, (13) 
^9 

with V being a space-time dependent element of the gauge group. In that case the gauge 
potentials even become pre-gauge invariant quantities. This class of solutions could also 
have been obtained by assuming the gauge fields to be spherically symmetric in the first 
place, which also leads to vanishing magnetic fields. However, as pointed out above, gauge 
potentials a priori have no meaning as observables and thus assuming them to be radi- 
ally symmetric is rather an ansatz than a necessary consequence imposed by a symmetry 
condition. 

With the particular form of the potential ([T3| ) it is now possible to find radially symmetric 
solutions of the Maxwell equations. For this purpose we note that the covariant derivative 
of a group valued quantity X can be rewritten in the form 

D,X = (d{VXV-^)^ V, (14) 
which suggests the introduction of the gauge equivalent quantities 

S = VEV-\ Q = Vpi^dV-\ W^ = VW^V-\ (15) 
Thus the Gauss-law and the zeroth component of (p^ 

[D, E] = pind, [v ■ D, = = -Ev (16) 
assume quasi-abelian form 

dS = Q, {v- d) W° = -Sv. (17) 

We may further choose the spatial components of v to point into radial direction, v = f/r. 
This entails that the spatial components of the induced current in ([TlD vanishes and the 
induced charge becomes proportional to W°, 

Q(r, t)=m^ j ^W°(r, t) = m^W\r, t). (18) 
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Due to the linearity of Eq. ([TTD and the particular form the diagonal color components 
of the electric field satisfy the differential equation 

(do + dr)(- + dr)S = m^S (19) 
r 

which will be in the center of our further investigations. 

We construct solutions in terms of eigenfunctions of the time derivative, S{r, t) = £{r)e^^* 
where the separation parameter uj may in general assume complex values. Eq. (|19|) can be 
straightforwardly transformed into the radial Schrodinger-equation for the hydrogen atom 

with K = a/ — ui"^ and p = nr. The solution of (|19D is composed of a singular and a 
regular part, 

E{r) = -—^ ^r(2 - i-)U{-l - -2, kt) + C2r^M{2 - 4, k, 



where U and M denote the confluent hypergeometric functions Ci,C2 are constants 
and the additional factor in front of the first term is chosen such that the series in r of the 
first summand in the brackets starts with unity. 

The quantity iuj/{K,p) which is the charge of the nucleus in the corresponding quantum 
mechanical problem, is an imaginary number for real u. This entails that k becomes purely 
imaginary when the first argument in the hypergeometric function M is a negative integer 
and its power series terminates. Consequently, no stationary regular bounded solutions - 
corresponding to the quantum mechanical bound states - exist. 

The singular part proportional to U which is usually excluded by the demand of regularity 
at the origin has in fact an interpretation as the electric field generated by a pointlike source. 
Let us 'measure' the charge corresponding to that field according to ([T7|). More precisely, we 
consider the divergence of S in the sense of distributions acting ClS db functional on C;;«(]R3) 
test functions ip, 

(Q, ip) = - [S, dip^ ■=- J dnj^ drr^S{r)dip{rn) (22) 

which by integration by parts leads to 

Vim{ATTr^£{r))ip{0) + J dil drr^ (^^-^f^) =■ <^) + (dS, v)- (23) 

The resulting charge is composed of a pointlike source cx 6^{r) which gets contributions only 
from the singular part of the electric field and the divergence which corresponds to the charge 
induced by the plasma. In what follows we will not consider the true vacuum solution oc M 
and concentrate on to physically more interesting situations with local charges. Particularly, 
we focus on an oscillating point charge and a charge which is constant for t < and vanishes 
for t > 0. 



6 



A. Monochromatic oscillations 



One particularly interesting class are stationary plasma waves induced by a periodically 
oscillating charge, Q}°'^^^{r,t;u)) = (5^(r)e*^* which electrical field is given by (|2T1) multiplied 
by e*^* with Ci = 1 and C2 = 0. Above the frequency uj = 2m, the inverse damping 
length K becomes an imaginary quantity, and we choose the root such that Q{k) < 0. The 
exponentially decreasing radial dependence in (^Tj) turns into a radially periodic one. The 
hypergeometric function U increases for large radial distances like r^+^^/'^ which combines 
with the prefactor to r~^^^^^'^ < r~^. 

For very large frequencies, /t —iu, the hypergeometric function independently from 
r becomes unity, and the radial dependence of the exponential vanishes. In this case, the 
electric fields oc exp{iu!t)/r'^ cease to propagate in the plasma. 

In the other extreme of a static charge, i.e. for a; — > 0, the electric field oc e~™''"(l+mr)/r^ 
corresponds to the exponentially screened Debye potential. 

It is instructive to consider the situation of perturbations with periods corresponding to 
mass threshold u. In that case the damping length 1/k becomes infinite and the argument 
in the gamma function in ( pT]) diverges. However, the infinities are just compensated by the 
hypergeometric function and one finds in the limit k —>■ 

f "=2™(r, t) = m' ^ e''^-^"ir3(2v^) (24) 

where x = mr, y = mt, and K3 denotes the Neumann function of imaginary argument. The 
large distance behavior is somehow in between the static exponential decrease and the high 
frequency power-law behavior. In fact, asymptotically expanding the Bessel function, 

^-=2™(r,t) = v^(2zx)5/^e2*^-"-2^ ( 1 + + o(- 

2x2 ' ' V 16v^ \x. 

^ ^i2^. f" exp (2.y - - 2(1 + .j^i - ^ + ^ . O (i)) (25) 

shows that the field is damped with an exponential oc ^/x only. The phase gets also contri- 
butions from ^/x, a constant term and inverse powers of the root. Asymptotically, constant 
phases propagate with velocity 

y ^^ = 2(1- x-'l^ + X-' + ^x-^'/^ + 0{x-^)\ 
dy \ 64 J 

which approaches twice the speed of light for x oo. Physical information about the 
propagation of waves in the plasma may be deduced from the acceleration of the phase 
velocity, 

^ = ~ 2x^3/2 - 6x-2 + ^x-^/^ + 0{x-'). 

(xy doc oz, 

The leading coefficient has a positive sign which means that a maximum of the oscillation 
runs slightly faster away from the origin than its neighboring maximum remoter from the 
source. A similar study in the regime x 1 shows that also there waves appear accelerated 
in positive radial direction with a maximal acceleration at x ~ 0.3 (Figl.). 
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B. Shock wave of a suddenly removed point charge 

In this section we study the situation where the locahzed charge is constant for t < 
and vanishes for t > 0. In a physical experiment, the charge may be a heavy quark which 
due to its large mass can be thought of as a well localized classical color source. At t = 
the quark may decay into a number of lighter particles which can escape the plasma or be 
thermalized depending on their energy, the temperature and the size of the plasma. In any 
case, we are faced with a scenario of a - with respect to the typical time scale m of the 
plasma - suddenly disappearing color-charge. 

To model the situation we have to find a superposition of the monochromatic oscillating 
solutions studied in the last section. This can be easily achieved by observing that the 
charge can be written as 

S\r)e{-t) = 5\t)^ [ du;^— = f / ^g^°-^(r, t; ^). (26) 

The path of integration is parallel to the real uj axes as shown in Fig. 2, and picks up 
the residue when closed in the lower half plane. The corresponding superposition of the 
monochromatic electric fields becomes 

= r^""- r(2-^)t/(-l-^,-2,/.r). (27) 

Prior to evaluating this integral by closing the contours above and below respectively, we 
have to define the root k = \/ 4m^ — a;^ for complex values of a;. It admits to branch points 
and it turns out to be most useful to chose the cuts parallel to the imaginary a;-axes. The 
sign is chosen that k — > ia; for large modulus of a; and |3fJ(a;)| > 2m. 

This choice admits to close the contour in the lower half plane for t — r < where the 
integrand vanishes exponentially on the arc. The only contribution comes from the residuum 
which gives the static limit 

^"''{r^) = 'Q^{\^mr) (28) 
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FIG. 2. Cut and pole structure in the w-plane of the decaying-particle solution. 



corresponding to the Debye potential. 

For t > r we may deform the contour to the arc oo, Q{uj) > on which the 

integrand vanishes but which picks up two cut-contributions. The contours around the 
cuts are well separated from the pole which facilitates considerably an asymptotical and 
numerical evaluation of the integration. In order to calculate the discontinuity we need the 
relation 

r(2 - a)U{-l - a, -2, z) - T{2 + a)U{~l + a, -2, -z)e^ 

aia^ — l)7r 



6 sin vra 



z-^M{2- a,A,z) (29) 



which can readily be derived from the integral representation of the hypergeometric function. 
Goosing the parametrization uj = 2m(— 1 + iu) for the left contribution, we get 



^f-yrf 2^3; /"°° , exp(-2iy + ix - 2yu + xu + xJu\Ju + 2%) 



127r 



smvr — 



\/u\/u+2i 
X — |— ui^ 

M{2 , 4, -2xv/m + 2i) + h.c. (30) 

V u + 2i 

where the complex conjugate part stems from the right cut. This expression is a particular 
superposition of regular solutions corresponding to vanishing local charge since the causal 
information of the decay propagates in the future directed light cone t > x oi the event at 
(x = 0,t = 0). 

The evaluation of the integral ( pOD is quite complicated since the integration variable 
appears not only in the argument but also in the index of the hypergeometric function. 



9 







FIG. 3. The electric shock wave propagating in the plasma. Below the light cone t < r the 
electric fields correspond to the Debye potential, are hit by the shock wave at t = r and vanish 
asymptotically for late times in the forward-directed light-cone. 

For values of x ~ 1 and larger, the situation is quite hopeless since even the asymptotic 
expansion of M for large arguments is spoiled by the diverging index at small u where the 
integrand contributes most. 

However, a numerical integration is still possible with the result shown in Fig. 3. The 
graph shows that waves are damped for large distances in direction orthogonal to the light 
cone. In the region close to the light cone r ~ i the amplitudes increase with the distance 
of the origin. 
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III. SUMMARY AND CONCLUSION 



We studied the spherically symmetric solutions of electric field configurations in a hot 
quark-gluon plasma. We showed assuming vanishing magnetic fields, which is a necessary 
compatibility condition with the symmetry imposed, that the electric fields which are a 
priori no physical observables in a non-abelian theory can nevertheless be given a physical 
interpretation. As in the abelian case, they arc radially directed vector fields. Their phys- 
ical components can be gauge transformed to depend on time and radius only. Using this 
property together with the equations of motions in the hard thermal loop approximation, 
we construct a differential equation for the fields. The solution of that equation consists of a 
singular and a regular part. Evaluating the corresponding charge in the sense of generalized 
functions shows that the singularity at the origin corresponds to pointlike charge, and the 
remaining part has an interpretation a induced charge. This admits to uniquely relate local- 
ized charge, induced charge and corresponding electric field in a manifestly gauge invariant 
way. 

We focus on two particular interesting physical situations: A periodically oscillating and 
a suddenly removed local charge. In the first case, we establish the electric field and concen- 
trate on oscillations with frequencies close to the frequency 2mDebye- The asymptotic field 
is found to oscillate with amplitude (mr)~^/^ e~^^/^, which is in between the exponentially 
damped static limit corresponding to the Debye potential and a power-law behavior for 
oscillations above cu — 2m. 

In the second example we study a suitably chosen superposition of the periodic solution 
which has constant local charge for t < and vanishing local charge for t > 0. The electric 
field is found to correspond to the static Debye case for t < r. In the causal region t > r, 
the electric field can be written as an integral which we evaluate numerically. The suddenly 
disappearing charge may be regarded as a model for a heavy, well localized quark placed 
in a hot quark-gluon plasma which spontaneously decays into a number of lighter particles 
which can escape the plasma or be thermalized. 

We found a gauge invariant concept to replace the QED-inspired idea of putting a local 
classical source into a plasma that potentially violates gauge invariance in the non-abelian 
case. We find that symmetry properties, and in particular spherical symmetry, already single 
out a class of solutions which admit a posteriori an interpretation of a localized charge. That 
result encourages to look for a symmetry-based definitions of physical observables, and in 
particular of the Debye mass. Work in that direction is currently under progress. 
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